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We characterize and construct time-independent Markovian dynamics that drive a finite- 
dimensional multipartite quantum system into a target (pure) entangled steady state, subject to 
physical locality constraints. In situations where the desired stabilization task can not be attained 
solely based on local dissipative means, we allow for local Hamiltonian control or, if the latter is 
not an option, we suitably restrict the set of admissible initial states. In both cases, we provide 
algorithms for constructing a master equation that achieves the intended objective and show how 
this can genuinely extend the manifold of stabilizable states. In particular, we present quasi-local 
control protocols for dissipatively engineering multipartite GHZ "cat" states and W states on n 
qubits. For GHZ states, we find that no scalable procedure exists for achieving stabilization from 
arbitrary initial states, whereas this is possible for a target W state by a suitable combination of a 
two-body Hamiltonian and dissipators. Interestingly, for both entanglement classes, we show that 
quasi-local stabilization may be scalably achieved conditional to initialization of the system in a 
large, appropriately chosen subspace. 
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I. INTRODUCTION 

Generating and manipulating highly non-classical en- 
tangled states in a robust scalable fashion is a central goal 
across Quantum Information Processing (QIP) [l| and 
quantum engineering [2| . While traditional circuit-based 
schemes are based on the application of suitable coher- 
ent control, implementing a desired sequence of quantum 
gates on a known fiducial state, progress in controlling 
open quantum systems is prompting a systematic explo- 
ration of incoherent control schemes using engineered dis- 
sipation. Access to dissipative control parameters, alone 
or in conjunction with Hamiltonian ones, may gener- 
ally allow for enhanced efficiency and added flexibility in 
achieving relevant QIP tasks. Notably, quantum simu- 
lation of arbitrary open-system dynamics on a quantum 
computer may be accomplished by suitable dissipative 
control on just one more "ancillary" qubit than needed 
to simulate closed-system dynamics l3| , as demonstrated 
in recent trapped- ion experiments [J, Q (see also [y]). 
Likewise, dissipative quantum-state preparation proto- 
cols may allow for arbitrary initial states to be driven 
into a desired target state by "all-to-one" control laws 
that have no counterpart for unitary evolutions, with the 



potential for inherent robustness against transient per- 
turbations. As a result, schemes for dissipative quantum 
control and generation of entangled states of light and 
matter via "reservoir engineering" have been extensively 
i_nvestigated theoretically in different contexts, see e.g. 
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As control methodologies that rely on engineered dis- 
sipation become increasingly attractive experimentally 
d, l26l - l28t . it is essential to develop rigorous control- 
theoretic characterizations of the set of target states 
achievable in open quantum systems given realistic con- 
trol resources, and to devise constructive control design 
protocols. In this context, an important task is provided 
by dissipative entanglement engineering in multipartite 
systems evolving according to a continuous-time quan- 
tum Markov process [2^. It is worth stressing that de- 
signing Markovian dynamics that admits a desired (pure 
or mixed) target state as its unique steady state (that 
is, achieves global asymptotic stability at the target) is 
relatively straightforward in the absence of control con- 
straints, and different constructive strategies exist. In 
fact, application of a single time- independent Hamilto- 
nian and Markovian noise channel suffices in principle in 
the generic case [ia| . Alternatively, if the system is com- 
pletely controllable in the absence of dissipation and the 
control actions can be enacted fast with respect to the 
noise time-scales, it clearly suffices to identify a purely 
dissipative dynamics that stabilizes any state with the 
same spectrum as the target, followed by rapid unitary 
control in order to suitably rotate such state into the 
target [ll|. 

In practice, however, physical evolutions arc not de- 



scribed by arbitrary Markovian master equations, and 
available coherent and incoherent controls are inevitably 
constrained. In particular, physically admissible Hamil- 
tonian and noise (Lindblad) operators are typically de- 
scribed by operators that act non-trivially on finite sub- 
sets of subsystems, that is, are Quasi-Local (QL) relative 
to the given tensor-product decomposition. Necessary 
and sufficient conditions for a pure entangled state to 
be stabilizablc in the absence of pre-existing "drift" dy- 
namics under QL constraints have been obtained in j36| , 
assuming purely dissipative means (so-called Dissipative 
QL Stabilization, DQLS for short). 

In this work, we continue and substantially expand 
our investigation, with the twofold goal of (i) obtaining 
a complete characterization of the DQLS setting in the 
presence of non-trivial drift dynamics [37|; (ii) defining 
and analyzing alternative QL settings, which may allow 
for a desired stabilization task to be achieved in cases 
where purely dissipative control is insufScient. While we 
shall in general allow for both time-independent dissi- 
pative and Hamiltonian (coherent) control resources, no 
access to a universal set of tunable Hamiltonians nor 
auxiliary quantum resources will be assumed, preventing 
switching schemes along the lines envisioned in [l8| or 
[24| from being viable. Specifically, two natural comple- 
mentary strategies will be pursued: to either augment the 
available control resources, by suitably combining (time- 
independent) QL Hamiltonian control with dissipation; 
or, if Hamiltonian control is not an option or is otherwise 
unfeasible, to insist on purely dissipative stabilization by 
restricting the initial condition to an appropriately cho- 
sen "attraction basin" . 

The content is organized as follows. We begin in Sec- 
tion III Al by introducing the relevant class of dynamical 
models and different QL stabilization settings, along with 
required known results on DQLS states. Hamiltonian- 
assisted and conditional QL stabilization are first ana- 
lyzed in Section IIIII and Section IIVI respectively, under 
the simplifying assumption that internal drift dynamics 
is absent or can be ignored. In both cases, we develop 
necessary conditions that the control Hamiltonian or, re- 
spectively, the initialization subspace must obey in or- 
der for the target state to be QL stabilizablc, and then 
present a constructive randomized algorithm that outputs 
a choice of stabilizing controls with unit probability. We 
reconsider and address the role of underlying Markovian 
drift dynamics in Section |Vl identifying necessary condi- 
tions for the desired stabilization to be compatible with 
the drift, and showing how different scenarios may then 
be tackled by adapting our previous results. Section IVII 
present explicit QL stabilization protocols for dissipative 
engineering of GHZ and W states in rt-qubit systems. 
While none of these highly entangled states can be gen- 
erated solely based on QL dissipative means, except for 
trivial exceptions [3^, we show that for both classes QL 
stabilization is achievable in a scalable fashion for a large 
subspace of initial states. For GHZ states, we addition- 
ally establish that no scalable procedure exists for achiev- 



ing stabilization from arbitrary initial states, whereas ar- 
bitrary initial states may be driven to a target W state 
by a suitable combination of two-body Hamiltonian and 
dissipators. Goncluding remarks are given in Section rVIII 
We include in separate Appendices the proof of conver- 
gence of our randomized construction as well as, for com- 
pleteness and direct reference, an algorithm for verifying 
global stability, originally developed in [35|, that is re- 
quired in the proof. 



II. PRELIMINARIES AND BACKGROUND 

A. Quasi-local quantum dynamical semigroups 

We consider a finite-dimensional multipartite open 
quantum system S consisting of n (distinguishable) sub- 
systems, defined on a tensor-product Hilbert space 

n 

v. = ^^V-a, a — 1, ... ,n, dimCHa) ~ da, dimCH) ^ d. 

a=l 

Let TiCH) denote the (convex) set of physical states (den- 
sity operators) of S, that is, trace-one, positive semi- 
definite operators with support on H. In analogy with 
probability distributions, we shall denote the support of 
p e ©("H) by supp(p) = range(p). It is then natural to 
use '3{H') to describe the set of density operators with 
support contained in a subspace %'(-%. 

In many QIP scenarios of interest, S undergoes 
continuous-time dissipative dynamics due to the coupling 
to a (physical or engineered) Markovian environment, in 
which case the resulting dynamics is described by a Lind- 
blad master equation [l|, |29| . In units where h ~ 1 . we 
may write 



= -^[H,p{t)] 



Y, (LkPml liLiLk 



(1) 
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in which C represents the most general form for the gen- 
erator of a completely-positive, trace-preserving Markov 
semigroup, {Tt = e^*,t > 0}, acting on 1){^). Here, 
H is an Hermitian operator associated with the Hamil- 
tonian of 5*, whereas the Lindblad (or noise) operators 
{Lfe} specify the non-Hamiltonian component of the gen- 
erator, resulting in non-unitary irreversible dynamics. It 
is worth recalling that, for a given Markovian generator, 
the decomposition into Hamiltonian and dissipative part 
in Ec^. ([l} is not unique, in the sense that 

C{H,{Lk}) = C{H,{Lk}), 

provided that the relevant operators are redefined as fol- 
lows, see e.g. Lemma 2 in |12| : 



Lk = Lk + CkI, H = H 



Y,{clLk-cuLl), CfeGC. 

k 

(2) 



Given a dynamical evolution of the above form, locality 
constraints are imposed by requiring both the Hamilto- 
nian and each Lindblad operator to have a non-trivial 
(non-identity) action only on certain subsets of subsys- 
tems, which may be distinguished by the geometry of the 
system and/or a physical coupling topology, and which 
we shall henceforth refer to as neighborhoods. Following 
[36| . neighborhoods {Afj} can be specified in full gener- 
ality as subsets of the set of indexes labeling the subsys- 
tems, that is, 

AT,- C {!,..., n}, j = l,...,M. 

Once a neighborhood structure is assigned onH, a list of 
reduced neighborhood states {p^f } can also be naturally 
associated to every state of S: 

p^^=TT^^ip), pe^in), j = l,...,M, (3) 

where Tr^. indicates the partial trace over the tensor 
complement of the neighborhood Afj, namely, Tij^. = 
®a(tj\f '^a- ^ Lindblad master equation as in Eq. ^ 
will be called QL if both its Hamiltonian and each of the 
noise operators are QL, according to the following: 

Definition 1 A Lindblad operator L is QL if there exists 
a neighborhood Afj such that: 

L = Lm^ ® /at, , 

where Lj\f. accounts for the action of L on the subsys- 
tems in Afj, and Ij^. := ^a4j\f ^a 'is the identity on the 
remaining ones. Similarly, a Hamiltonian H is QL if it 
admits a decomposition into a sum of QL terms: 



It is immediate to verify that the QL property is well 
defined with respect to the freedom in the representation 
of the Markov generator, since the transformation in Eq. 
^ preserves the quasi- locality character of H and {Lk]. 
We also note that introducing locality constraints based 
on neighborhoods formally encompasses different specific 
notions that are encountered in the physical and QIP lit- 
erature [ll [nl, [ij] , where locality is typically associated 
with operators of fixed maximum weight t (also called "i- 
body" or "t-local" interactions) and/or "distance" on a 
graph. For instance, when neighborhoods coincide with 
individual subsystems, strictly local (or "single-site") dy- 
namics is enforced, with the corresponding propagator 7t 
being fully factorized relative to the multipartite struc- 
ture, Tt = ®a=iT'a,t ll2l- Likewise, for i = 2, allowing 
for arbitrary two-body (possibly long-range) interactions 
correspond to identifying neighborhoods with all possible 
subsystem pairs, whereas interactions between nearest- 
neighbors sites may be accounted for by restricting the 
neighborhoods accordingly. 



B. Stability definitions and standard form for 
stabilization 



As mentioned, our focus in this paper is on stabilizing a 
desired pure state of S. The relevant notions of stability 
are given in the following: 

Definition 2 Let pd = |*)(*| e S)('H) be a pure state of 
S evolving under Eq. ^ . 

(i) Pd is Globally Asymptotically Stable (GAS) if for 
every initial condition po € 'S{'H) we have 



lim e^* 



[Po] = 



(4) 



(ii) IfT-L'^T-L is a proper subspace ofH, pd is Condi- 
tionally Asymptotically Stable relative to T-C (or simply 
n'-AS) if Eq. 0) holds for pa e '!){%'). 

A necessary condition for a state to be GAS is that 
it is invariant |5l| . namely, that it is in the kernel of 
the Liouvillian, C{pd) = 0. The following proposition, 
that follows from the general results of [IJl, provides a 
particularly simple way to check for invariance if pd is 
pure: 

Proposition 1 Let the dynamics be driven by 
C{H, {Lk}) as in Eq. |ip. Then a state 

Pd = |'J')(^| G 1){'H) is invariant if and only if 



Lfe|^)=4|*), 4eC, VA:, 



(5) 



\Xd 


Xp] 


[Xq 


Xr\ 



H\^) = h\^), H = H-^ llMLk ~ IkLl), h e m 



Proof: Let Hd = span{|^)} and, with respect to the 
natural representation induced by the partition H — 
Hd ffi Hd , let operators be given in the block form 



X 



where the top- left block is one-dimensional. Given Corol- 
lary 1 in [l2|, it follows that pd is invariant if and only 
if ifcl*) = lk,d\'^) = lk\^) for aU fc, and the entries in 
the F-block olDcy iHp - (1/2) J^k ^kLp,k = 0. Clearly, if 
£fe = in Eq. ©, then it follows that Hp = = Hq hy 
Hermiticity. Therefore, H = H is block-diagonal, with 
H\-^) = hd\-^) = hd\-^), as stated in ®. If 4 7^ for 
some fc, we can use Eq. ([2]) to redefine Lk = Lk — tkL 
and, accordingly, H as in ©. Since now £k = 0, H is, 
again, block-diagonal, with |^) being an eigenstate with 
eigenvalue hd = h, as stated. D 

Thanks to the above result, it is always possible to ex- 
press a Markovian generator that asymptotically stabi- 
lizes a pure state in a standard form, in which the desired 
target state |^) is annihilated by all the noise operators 
(thus being a "dark state" in quantum-optics language 
[III) and is an eigenstate of the corresponding Hamilto- 
nian. That is: 



Corollary 1 If a generator C{H, {ife}) makes pd ~ 
|\I/)(^| G 'D{'H) GAS, then the same generator can be 
represented in a standard form C{H, {Lk}), in such a 
way that H\^) = h\^) and Lfcl*) = 0, for all k. 



C. Quasi-local stabilization settings 

In general, a QL stabilization problem will entail the 
specification of the desired target state and the relevant 
QL constraints, along with a description of the target dy- 
namical model and control parameters. While, as men- 
tioned, we shall focus here on pure target states, different 
stabilization settings may be envisioned depending on the 
nature of the available control resources and initialization 
capabilities, and on the existence of non-trivial Marko- 
vian dynamics in the absence of control, hereby referred 
to as drift in the standard control-theory terminology and 
characterized by a Liouvillian generator Cq . Specifically, 
let us consider a class of dynamical models as in Eq. ([T]) , 
with the corresponding generator C = C{H, {Lk}) taken 
to be in standard form henceforth, and 



p{t)^^p(t))^[Co + C,]{pit)). 



(7) 



Here, £o{Hq,{L^}) and Cc{Hc,{Dk}) are associated to 
the drift and control dynamics, respectively, and H = 
Hq + Hc, {Lk} = {ifci L>k}- Hamiltonian control is intro- 
duced by the time- independent He, whereas each of the 
Lindblad operators Dk represent a different incoherent 
control knob on S. We can then introduce the following 
definitions. 

Defimtion 3 A pure state pd = |*)(*| e ©("H) is: 

(i) Quasi-Locally Stabilizable (QLS) if there exist a QL 
Hamiltonian He and QL noise operators {Dk}k=i....,K 'in 
standard form, such that pd is GAS for the controlled 
Lindblad evolution in Eq. (f^. 

(a) Dissipatively Quasi-Locally Stabilizable (DQLS) 
if it is QLS with He = and QL noise operators 
{Dk}k=i,....K in standard form. 

(Hi) Conditionally DQLS with respect to a subspace 
W C. % (or simply "H'-DQLSJ, if there exist QL noise 
operators {Dk}k=i,....K in standard form, such that pd is 
conditionally AS relative to %' for the dissipatively con- 
trolled Lindblad evolution in Eq. (f^. 

Two features follow from the global stability properties 
given in Definition [21 while exact QL preparation cannot 
be achieved in finite time, convergence happens exponen- 
tially fast in time |35| . Furthermore, if the target state is 
stabilizable, according to the appropriate definition, the 
stabilizing dynamics is robust with respect to the initial 
state [12j thus, equivalently, the control parameters are 
all-to-one [18|. 

Clearly, the QLS setting includes both purely dissi- 
pative and conditional stabilization as special instances 
where either (or both) the admissible control actions and 
initial states are restricted, respectively. Introducing a 



separate DQLS notion is motivated by the fact that it 
allows for a simpler mathematical treatment (see Sec. 
Ill Dp while being adequate for important classes of en- 
tangled pure states and, from a practical standpoint, 
potentially simpler to implement. For states which are 
not DQLS, QLS may or may not be achievable. Even 
in the latter case, conditional stabilization may offer a 
practical method of choice in situations where a two-step 
"switched" dynamics, that first initializes 5* to the in- 
tended subspace, and from there enacts dissipative stabi- 
lization, is preferable over identifying and implementing 
a combined Hamiltonian-dissipative control action - for 
instance thanks to the presence of a conserved quantity 
[5[ . While a number of illustrative examples will be dis- 
cussed in Sec. I VII ensuring that the required number K 
of noise operators scales favorably with the complexity 
of the neighborhood structure is crucial in practice. The 
situation is straightforward for the DQLS setting, since it 
follows directly from Theorem [1] below that a single noise 
operator per neighborhood always suffices [K = M). In 
general, the randomized algorithms that we will provide 
to achieve stabilization in case (i) and (ii) will also out- 
put, when successful, a noise operator per neighborhood, 
as we shall see. 

We begin by presenting two preliminary results that 
directly generalize their counterparts in Ref. [36| for the 
DQLS case. First, let U = (3)"=! Ua be an arbitrary 
local unitary (LU) transformation. Then it is straight- 
forward to see that the above stabilization notions are 
invariant under arbitrary LU transformations of the tar- 
get state, as desirable given that entanglement proper- 
ties are themselves preserved under LUs. The following 
Proposition may be established through the same steps 
used in Lemma 2.2 and Proposition 2.3 in p36| : 

Proposition 2 If pd is QLS and U is any LU, then p'^ — 
U pdU'^ is also QLS. If Pd is T-L'-DQLS and U is any LU 
that leaves H' invariant, then p'^ = U pdW^ is also T-L' - 
DQLS. 

A second basic yet useful result regards the structure 
of the operators {Dk}: the Lemma that follows shows 
that the support of a QLS state must still be, as in the 
DQLS case, contained in the kernel of the noise operators 
written in standard form: 

Lemma 1 Assume that the generator associated to QL 
{H,Dk} leaves pd ~ |'I')(^| invariant. Then, for each fc, 
we have supp(pA/k) ^ ker(DA^j^). 



Proof: If the state is invariant, by Proposition [I] |^) 
must be in the kernel of each Dk. Thus, with respect to 
the decomposition H = Hd ffi "H j" , with Hd = span{ | ^I*) } 
as before, every Dk must be of block form |13j : 



Dk = 



Dp,k 
DR,k 



which immediately imphcs DkPdD\ ~ 0. It then fol- 
lows that Trjy-j^ {D^pdO].) = 0, therefore Tr_;\/-^ {Dj^^ <g) 
Ij^^PdDj^ (g) Ij^^) = 0. Thus, it also follows that 
Dj\ff.pj\f^Dj^ ~ 0. If we consider the spectral decomposi- 
tion pj\fj, = 22i Qj\'f'j){4'j\j with Qj > 0, the latter implies 
that, for each j, D^i^\(j)j) {(j)j\Dj^ — 0. Thus, it must be 
supp(pA^^) C ker(D7v'k)i as stated. D 

For added clarity and notational simplicity, we shall 
from now on assume that the dynamics of the system 
is drift-less, that is, £o = in Eq. ([7]). After briefly 
recalling the main results on DQLS, we will present our 
new results on Hamiltonian-assisted and conditional QL 
stabilization under the simplyfing drift-less assumption, 
in Sec. IIIII and Sec. IIVI respectively. We will then allow 
for drift dynamics and explicitly address its role in Sec. 

m 



D. Prior results for quasi-local dissipative 
stabilization 

A characterization of DQLS states, leading to a simple 
linear- algebraic algorithm to test whether a given pure 
state is DQLS, may be obtained based on the proper- 
ties of the reduced states on the neighborhoods [52|, as 
defined in Eq. ([3]). Let 

m = supp(p^, ® /^j, no = f]ni. (8) 

k 

We thus have the following: 



state of a frustration-free parent Hamiltonian. Then the 
QL structure of H may be naturally used to derive a 
stabilizing semigroup [11|. It is easy to show that this 
condition is also necessary, leading to the following: 



Theorem 1 |36| A pure state pd = |*)(*| is DQLS if 
and only if 



supp{pd) = Ho- 



(9) 



The proof of this Theorem includes an explicit con- 
struction of a choice of stabilizing noise operators, with 
their general block-structure highlighted in Lemma [T] 
It also points to natural connections with the formal- 
ism of parent Hamiltonians and the concept of frustra- 
tion from many-body physics. Let H = '^j.Hk, with 
Hk = Hj^^ (X) Ijij- and, as before, the index k refers to the 
kth neighborhood. A QL Hamiltonian H is said to be a 
parent Hamiltonian for | ^P) , if | ^) is an exact ground state 
for H [30|. In addition, |^) is said to be frustration-free 
if it is the exact ground state of each Hamiltonian Hk 
separately, that is, {^\Hk\'^) = min \{Hk), Vfc, where 
A(-) denotes the spectrum of a matrix. Remarkably, it 
is known that for a large class of product entangled-pair 
states (PEPSs) or, in one spatial dimension, inatrix prod- 
uct states (MPSs), a frustration- free parent Hamiltonian 
may be constructed [5311 which has the desired state as its 
unique ground state J33 - |34| . The QL properties of this 
Hamiltonian are determined by the (minimum) "bond di- 
mension" of the corresponding PEPS or MPS represen- 
tation. Suppose that a pure state is the unique ground 



Corollary 2 |3i| ^ state pd = |*)(*| is DQLS if and 
only if it is the unique ground state of a frustration-free 
parent Hamiltonian. 

Despite the formal points of contact, it is crucial to 
remark that in MPS-based stabilization approaches [il| 
the relevant QL notion is intrinsically state- dependent, 
whereas it is taken to be a fixed problem input in our 
control-motivated approach. From a practical stand- 
point, the DQLS class includes important representative 
sets of entangled states - most notably, all stabilizer and 
graph states [111 |36|, relative to the natural choice of 
neighborhoods associated with connected nodes on the 
graph. Still, paradigmatic examples of genuinely multi- 
partite entangled states such as GHZ and W states can 
be easily seen to fail the DQLS test in Eq. (jH]) , except in 
"fully connected" (for example, "star") coupling topolo- 
gies where the QL constraint becomes effectively trivial. 
While we defer to [3g and Sec. IVII for further discussion, 
we shall proceed to separately formalize and analyze the 
general QLS and the conditional DQLS scenarios next. 



III. QUASI-LOCAL STABILIZATION WITH 

CONSTANT DISSIPATIVE AND HAMILTONIAN 

CONTROL 

A. Necessary conditions 

Assume that the target state pd = |^) (*| is not DQLS, 
and let Hd = span{|\E')} as before, corresponding to the 
support of the desired state. In terms of 'Ho, defined as 
in Eq. ([8|), this translates to: 

dim (-Ho) = dim(p|-H^.) =do>2. 

k 

Our first result is the following necessary condition for 
QLS (recall that H = He in the drift-less scenario we 
consider for now): 

Proposition 3 pd is QLS but not DQLS only ifHd is an 
invariant subspace for the Hamiltonian He and no other 
invariant subspace is contained in (or equal to) H-q. In 
particular, one can choose H^ so that i?c|^) = 0. 



Proof: If there were another invariant subspace for the 
Hamiltonian with support in Hq, the latter would be, 
by definition of Hq, also in the kernel of each Dk, and 
hence it would be invariant. Invariant subspaces always 
contain at least an invariant state. It then follows that 
Pd could not be GAS. If i^d*) = A|*), A 7^ 0, we can 



always choose H'^ = He — XI instead, which is also QL if 
He was. D 

When such a stabilizing Hamiltonian He exists, one 
must look for noise operators {Dk} such that "Hd is 
the only invariant subspace for the whole generator 
£c{Hc,{Dk}), and thus makes p^j GAS. Lemma [1] sug- 
gests that the most effective choice of noise operators 
can stabilize "Ho, but do no better than that. In order 
to specify what the action of a stabilizing Hamiltonian 
would be, it is convenient to pick an orthonormal basis 
for Ha, which includes the target state: 



Ho = span{|*),|$i) 



I*.)}, 



= rfn - 1. 



One would hope that Hc\^j) ^ Ho for each j. However, 
fulfilling these conditions is clearly not necessary, and 
in fact it need not be possible given the QL constraint. 
However, in the simplest case of do = 2, the above idea 
leads to a specialized formulation of Proposition [3] 

Corollary 3 pd is QLS but not DQLS, with Ho ~ 
spanjl^), |$i)}, only if there exists a QL Hamiltonian 
Hr such that 



Hem = 0, 



i^c|*i> i Ho 



(10) 



Recall that in the DQLS context (Corollary [2]), a QL 
parent Hamiltonian was naturally associated to the ac- 
tion of the noise operators, in such a way that Ho was 
the common ground eigenspace of all the QL components, 
that is, with no frustration involved. In contrast, the nec- 
essary conditions provided above clearly show that an ef- 
fective control Hamiltonian He cannot he frustration- free 
with respect to Ho, since it must destabilize some part 
of Ho in order to attain GAS of the target state. 

In Sec. IVII we will employ Corollary [3] to construct 
a stabilizing Hamiltonian for both GHZ and W states 
which, as noted earlier, are newer DQLS under non-trivial 
QL constraints. 



B. Randomized construction of stabilizing 
Hamiltonian and dissipators 

Even if we succeed in finding a Hamiltonian He that 
satisfies Proposition [3] (or Corollary [3]) , a procedure for 
determining the existence and, possibly, the actual form 
of the stabilizing Dk is required in order to establish QLS. 
As we will illustrate with an example in Sec. IVIBl not 
all choices of {Dk} satisfying Lemma{l\are effective. In 
fact, care is needed in ensuring that the interplay between 
Hamiltonian and dissipative control introduces enough 
"mixing" and does not allow for other invariant sets to 
exist. 

With that in mind, we shall invoke a randomized design 
ap pro ach that expands our earlier use of randomization 
in [35| to both the Hamiltonian and the dissipative com- 
ponent, and prove that a generic choice of noise opera- 
tors which satisfy Lemma [T] and stabilize Ho will suffice 



to achieve QLS. Specifically, assume that we represent 
the desired stabilizing operators He = X]fc ^k, and {Dk} 
in parametric form: 



^^■ = E' 



^jkCTjk, 



D, 



= E/^. 



jkO'jk, 



(11) 



where ajk,Pjk S K are chosen at random with uniform 
distribution in an interval X = [—7,7] of the real axis, 
and {ajk} is a basis [5J] for the space of QL operators 
on the fcth neighborhood A/fc- Our main result is then 
contained in the following: 

Theorem 2 // there exists a choice of ajk, /3jk G I that 
makes pd QLS, then almost any choice of Ujk, jSjk G I, 
that makes pd invariant, makes it QLS as well. 

While the proof of the above Theorem is rather tech- 
nical (see Appendix I VIip . the meaning is clear: imposing 
invariance requires "fine-tuning" of the parameters, how- 
ever once invariance is ensured, if GAS is possible, then 
it comes almost always for free. An algorithm for con- 
structing a stabilizing QL Hamiltonian and achieve QLS 
may then be provided as follows. 

Step 1: Imposing quasi-locality of {He}. Pick a 
product operator basis for the d^-dimensional 
space of linear operators 23 (H), say, {at^ (g) . . . (g) 
cri^\ia = l,...,dl}. Let He = Y^k^k, and as- 
sociate each Hk and Dk to a vector hk and dk, 
respectively. Let D = [di\---\dn] be the matrix 
of the coefficient of the noise operators relative to 
the above product basis, and Bk the orthogonal 
projection onto the subspace generated by the ba- 
sis elements that are QL with respect to Mk , with 
Bj: = I — Bk- We must then require, for each k: 



Bihk=i), Bid. 



kO-k 



0. 



Step 2: Ensuring invariance of Ho- Impose the linear 
constraints on the QL noise operators and on the 
Hamiltonian, namely Dk\ip) =0, i?c|V') =0 ^oi all 
k. If Pq is the matrix representation of the linear 
(super)-operator Po{X) = A|-0)(V'| with respect to 
the the chosen basis, this is equivalent to require 
that for each k. 



Poh = 0, Podk = 0. 

These two steps translate, for each k, in the follow- 
ing homogeneous systems of linear equations: 



B 



k 

Po 



dk = 0, 



B 



k 

Po 



hk = 0. 



(12) 



Step 3: Enforcing convergence by randomization. If 
the constraints above allow for non-trivial solutions 
(and hence a subspace of solutions since the system 
is homogenous), choose the free variables uniformly 
at random within a finite interval [—7, 7] C M. 



We thus have the following immediate corollary of The- 
orem [51 

Corollary 4 // pd is QLS, the QL generator 
Cc{Hc, {Dk}) constructed in Steps 1-3 makes it 
GAS with probability one. If pd is not QLS, then it is 
not stabilized by the constructed Hamiltonian He and 
noise operators {Dk\. 

Since checking whether a state is the unique equilib- 
rium of some Lindblad dynamics is straightforward {e.g., 
by checking that the superoperator form of the corre- 
sponding generator has a unique unit eigenvalue), the 
above algorithm can be used as an explicit test for QLS. 



IV. CONDITIONAL QUASI-LOCAL 
DISSIPATIVE STABILIZATION 

A. Characterization of conditional stability 

When a state is not DQLS and Hamiltonian control 
is not viable, we can return to a purely dissipative con- 
trol setting and analyze our second proposed stabilization 
strategy, namely to restrict the initial state to a given ini- 
tial subspace of T-L. With T-Ld = span{|^)} and T-Lo given 
in Eq. ([5]) as before, let us additionally define the follow- 
ing subspaces: 



Hence p € 'D{H') would not converge to pd- 



a 



Hw — T~Lq G T-Ld, 



rir — ri G riQ- 



By construction, Hw corresponds to the portion of the 
subspace stabilizable by purely dissipative means that is 
orthogonal to the target. Assume that the dynamics is 
given by Eq. (O, with £ = 0, associated to QL noise 
operators {Dk} that make S)('Ho) GAS. We aim to char- 
acterize which choice(s) of H' can make pd conditionally 
DQLS (with the obvious requirement that H' 2 Hd)- 
A first necessary condition is provided by the following 
Lemma: 

Lemma 2 pd can be T-L' -DQLS only if 

W CHeHu,=Hd<S)[HeUa]. (13) 



Proof: Since 1){T-Lw) £ ©("Ho), and by hypothesis 
'D{'Ho) is invariant under the Lindblad dynamics Cd in- 
duced by the {Dk}, then 'D{Hw) is invariant as well. 
If "H' C H G Hw did not hold, then there would be a 
p e 2)(?^') such that Tr {pU-uJ > 0, with H«^, denot- 
ing the projector onto Hw- Define Pw = Tr (H-^^p) > 0, 
Pw = p~'^Il-H^, pU-H^ e "^iHw), and write p = p^Pw + ^P, 
for some Hermitian but not necessarily positive opera- 
tor Ap. Since the dynamics is linear and pw invariant, 
Tr(pH-H„) cannot decrease along the dynamical flow |13 |. 
therefore we may write 

lim Tr(e^"*[p]H„ J = Tr(e^"*[Ap])+p„Tr(p^) > p„ > 



The following "enlargement lemma" indicates that we 
can construct larger stabilizing subspaces out of smaller 
ones, both weakening our constraints and gaining a key 
property: invariance. Formally: 

Lemma 3 // pd is W'-DQLS, it is also W -DQLS for 
some %' such that %" C %' and 1)(T-L') is an invariant 
set. 



Proof: In order to make S)(H') invariant, it suffices 
to combine the supports of all trajectories originated in 
1){J-L"). Given two subspaces Hi, 2, let us denote by 
T-Li V H2, the smallest subspace that contains both of 
them. Define 



T-Ln 



\/supp(c^°'[po]), W ^ V Up 



PoGSCH") 



By assumption, Ud C U' , U" C %' and ©("H') is in- 
variant. Lastly, all trajectories originated in Ti{T-L') can 
be obtained from trajectories originated in 1){T-L"), which 
by hypothesis all converge to pd, establishing the desired 
result. D 

Motivated by the above, we now restrict to subspaces 
T-L' such that 'D{T-L') is invariant. With a slight abuse in 
terminology, we shall simply refer to subspaces obeying 
such a property as invariant. We are then ready to state 
our main result: 

Theorem 3 Assume that the QL dissipative dynamics 
generated by the {Dk} makes ©("Hq) GAS. IfTi' satisfies 
the necessary condition given in Eq. (|13p and is invariant 
for all {Dk}, then pd is Ti' -DQLS. 

Proof: Since Ti' is invariant for all {Dk}, it follows that 
'p(T-L') is a positive invariant set. By LaSalle theorem 
[3l| . this means that all the trajectories starting in ©(W) 
converge to its largest invariant set. On the other hand, 
by hypothesis, S)('Ho) is GAS, so the largest invariant set 
must be contained in 2) (Ho) H T){Ti'). Since by p^ we 
know that Tiw -L T-L', it follows that 



S)(Ho)nD(H') = s(Hd) = Pd, 



as desired. 



D 



A further advantage of considering an invariant sub- 
space Ti' is highlighted in the following corollary: if ini- 
tialization in Ti' is faulty, the error on the asymptotic 
result remains upper-bounded by the preparation error: 

Corollary 5 Let pd be Ti' -DQLS, with Ti' invariant un- 
der the dissipative dynamics. If p is such that 



then we have: 



1 - Tr(n«,p) = e > 0, 



lim Tr(e^°'[p]n«0 < e. 

i— >-oo 



Proof: Since H' is invariant, Tr(n-H'p) is non decreasing 
along the trajectories e'-'^p, for t > 0, see [12|- Q 

While it seems hard to devise a fully general strategy 
for finding good choices of %' and associated QL dissipa- 
tors so that pd is "H'-DQLS, constructive results may be 
obtained if the problem is further constrained. In par- 
ticular, note that from a stabilization point of view, two 
subspace decompositions play a specially important role: 

(i) one associated with the initial (i = 0) state-space 
structure, V. = %' (B W-^, with H' containing the states 
to be attracted toward pd', 

(ii) one associated with the final (t — )> oo) state-space 
structure, Hq ~ Hd ffi "H^,, with Hw containing the un- 
wanted states toward which dissipative stabilization oc- 
curs. 

In a way, the two classes of conditional stabilization prob- 
lems we solve can be seen to arise by imposing some 
natural constraints on H'-^ and Hw, in case (i) and (ii), 
respectively. Either way, it is worth remarking that any 
state in '£){%') is asymptotically converging to pd'- if con- 
trol capabilities are enough to prepare a subset or even 
a single state in this set, dissipative QL preparation can 
be achieved. We begin by addressing case (ii), which is 
directly motivated by W states and is technically simpler. 



Conditional stabilization under constraints on 
the final attractive set 



Specifically, our candidate subspace for conditional in- 
variance is defined as follows: 



By Lemma [T] and the definition of T-Cj^ , T-Lw Q C\k'^k- 
Hence, %' obeys the necessary condition established in 
Eq. dni: 



W CHeU^ 



(15) 



We now need to construct noise operators Dk ~ Dj^^ (E) 
Ij^^, such that "Ho and H' are invariant. To this aim, 
define HV, = "H^^ G H^^ , which by Eq. 1^ is not 
empty, and consider Dj^^, with the following block struc- 
ture, with respect to the decomposition of the QL space 



Da/, 




Dp,A^, 
i?flx. 



(16) 



The above structure can be exploited, e.g. via a choice 
of Dp^^^ , Dfi^_\f^ of "ladder form" [l2| , to render each of 
the H]^ (hence Ho) GAS. We hereby consider this or 
an equivalent choice, so that Hq is made GAS with QL 
operators. In addition to this, the block structure in Eq. 
(|16p ensures the required invariance property of H': 

Proposition 4 A choice of Dk satisfying Eq. (fT6|) en- 
sures that H' — Hq PI^ Hk is invariant. 



In order to formulate and interpret the required prop- 
erty that Hw must obey to allow for H' and associated 
QL dissipators to be systematically constructed, some 
additional definitions are needed. Let 

'^If, = supp(pA^J, H'xf^ = H^/^ © Hlf^, 



?^^^ = supp (Tr^jn^ J) , 

where as before Hn^ is the orthogonal projector on Hw 
By construction, Hw C Hq and H"^^ C H]^^. As we 
shall establish, the key property that Hw must obey is 
the following strict inclusion for each neighborhood: 



n^, cH 



■A/'fc- 



(14) 



Some intuition on the above requirement may be built 
as follows. Wc know that Hj\f^ must be in the kernel of 
all the noise operators D^, and hence the dynamics in 
this subspace is trivial. We thus would like to be able 
to distinguish Hd from Hw, and "push" the evolution 
towards the former but not the latter. If Eq. (fH|) is 
obeyed, we know that if a pure state is in Hq but has 
no QL support on any of the HJj- , then it must be in 
Hd and hence in pd. Formalizing this intution provides 
us with a way to construct a suitable subspace H' and 
utilize Theorem [3] to prove convergence. 



Proof: Given the matrix structure, we have: 

DkH' QDkiHBHk) 

= (D^. ® I^, ) {mk ® ^k ) ® ^AT, ) (17) 
c(hV,©h;,J©H^,. 

Furthermore, 
f]HkCHk^ H-^^QHj^^ ± (Hk®^k)®^Ar,- (18) 

Thus, by combining Eqs. p7p and p^ . we obtain 
DkH' -L Pl^ Hk, so by definition of H' wc have DkH' C 
H'. D 



Given Eq. (|15p and the above proposition, we can 
then apply Theorem [3] to establish that the constructed 
dynamics makes pd conditional T^'-DQLS. We summarize 
the results of this section in the following: 

Corollary 6 Let pd and the given neighborhood structure 
{A/fe} fee such that Eq. p4p holds. Then by choosing 
H' ^He nfe('^k ® ■^atJ; and a set of QL {Dk] that 
satisfy Eq. ^, pdis H'-DQLS. 

As remarked, a notable example of states satisfying the 
property in Eq. (fTi|) is the class of W states. Their stabi- 
lization using this technique will be explicitly addressed 
in Sec. IVll 



C. Conditional stabilization under constraints on 
the initial attraction basin 

In this case, the additional assumption we impose on 
the noise operators {Dk} is that both %' and H'^ are left 
invariant. This enables us to test a candidate subspace 
H' and construct associated QL {Dk} by employing a 
variation of the randomized algorithm presented in Sec. 
IIIIBI for general QL stabilization. 

Technically, the key simplification that the above in- 
variancc requirements translates into is the fact that only 
linear constraints are imposed on the noise operators. 
This may be seen by writing the block decomposition of 
the Dk with respect io'H = 'H' ® V.'^: 



Dk = 



Ds.k Dp,k 
Dqm Dr^u 



In order for %' to be invariant, it must hold that nj 



E^ 



DlkDpM 



(19) 



and Dq^k = for all k. In terms of the coefficients jSjk 
that parametrize Dk as in Eq. (fTTj) . the above Eq. (fT9|) 
is clearly a set of quadratic equations. However, if T-L'^ is 



required to be invariant as well, then Eq. p^ becomes 
Dp^k = and all the constraints are linear. Hence, the 
proof of Theorem [5] (see Appendix IVII[) carries over to 
this case upon restricting to the subspace of fijk such 
that pd, H' , and "H'-^ are invariant. Let ajk = for all 
j,k. We thus have the following: 

Theorem 4 // there exists a choice of /3jk that makes 
Pd H' -DQLS and H invariant, then almost any choice 
such that Pd and Ti,'^ are invariant makes pd T-L' -DQLS 
as well. 

We can thus provide the following randomized condi- 
tional stabilization algorithm: 

Step 0: Finding a candidate %' . Given Lcnimi [2] and |3l 
we need T-L' D T-Ld such that T-L' -L Tiw. The largest, 
obvious candidate is of course Ti' ~ TiQTiw A trial 
choice may otherwise be dictated by physical con- 
siderations and available experimental capabilities: 
e.g., in situations where a conserved quantity asso- 
ciated to an observable O exists, such that Tid and 
Tiw belong to two different eigenvalues of O, a way 
to construct "H' is to identify it with the eigcnspace 
containing Tid. 

Step 1: Imposing quasi-locality of {Dk}. Pick, as 
before, a product operator basis for ^{Ti), say, 
{o-j^ (8) . . . (8) cri„ \ia ~ 1, ■ • ■ , d^}. Let Dk be associ- 
ated to a vector dk, whose components represent 
our free design parameters. Impose the QL re- 
quirement by demanding that the components as- 
sociated to basis elements that are not QL van- 
ish: let Bk be the orthogonal projection onto the 



subspace generated by the basis elements that are 
QL with respect to A/fc , and B^ = I — Bk- Write 
D = [di\ ■ ■ ■ \dn] for the matrix of all the free pa- 
rameters in compact form. We must then require, 
for each k: 

Bidk = 0. 

Step 2: Ensuring invariance of Tio- Impose the linear 
constraints Dk\ip) ~ for all k. If Pq is the ma- 
trix representation of the linear (supcr)-opcrator 
Pq{X) = X|?/»)(-i/'| with respect to the the chosen 
operator basis, this is equivalent to require, for each 
k: 

Podk = 0. 

Step 3: Ensuring invariance of Ti' and Ti'^. We now 
impose a block-diagonal form for the Dk with re- 
spect to the decompositon Ti =T-L' ®T-L'^ . That is, 
upon writing 



Dk = 



Ds.k 
Do.k 



Dp,k 
Drm 



we require that |12l |: 

Dn.k^Q, Dp.k = 0, 



(20) 



for all k. Note that Eq. (pil)) is a linear constraint, 
and can be imposed by requesting: [Z?fc,H'] = 0, 
where H' is the orthogonal projector associated to 
Ti' . In terms of the vectorization employed in the 
previous steps, this can in turn be rewritten as: 



P'dk = 0, 



(21) 



with P' denoting the matrix representation of the 
super-operator associated with the commutator 
with H', P'{X) = [X,U']. For each k, the above 
(three) steps translate in the following (homoge- 
neous) systems of linear equations: 



Cdk 



B 



k 
P' 



dk^O 



(22) 



If a non-zero solution exists, then a subspace of 
solutions dk exist, each corresponding to a QL op- 
erator Dk. 

Step 4: Forcing convergence. Choose the free vari- 
ables uniformly at random within a finite interval 

[-7,7] CR. 

We thus have the following Corollary as a direct con- 
sequence of Theorem 21 
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Corollary 7 Let pd denote the target pure state. As- 
sume that there exist an invariant subspace %' and QL 
{Dk} such that pd is H'-DQLS and T-L'^ is invariant. 
Then the above Steps 0-4 fii^d a choice of %' and {-Dfc} 
with probability one. 

From a practical standpoint, Proposition[7]providcs an 
explicit algorithm to test whether a candidate subspace 
%' makes the target pd H'-DQES and, if so, it also out- 
puts a choice of effective noise operators {D^}. 



QUASI-LOCAL STABILIZATION WITH 
DRIFT DYNAMICS 



In realistic scenarios, the uncontrolled system may 
evolve under an internal QL dynamics, due to a pre- 
existing Hamiltonians and/or couplings with Markovian 
environments, resulting in a non-trivial, known drift gen- 
erator, Co{Ho, {L1}) ^ 0, in the notation of Section lllCI 
We first discuss the simpler case where only Hamiltonian 
drift is present (L° = for all k). 



A. Drift Hamiltonian 

We preliminarily note that if complete QL Hamiltonian 
control is available over S with respect to the same local- 
ity notion of Hq or a less restrictive one, then a control 
Hamiltonian He = —Hq can be applied to undo the ac- 
tion of the drift, given that Hq is known. In this case, 
the problem becomes again effectively drift-less, and the 
results developed in the previous sections apply. We thus 
assume here that limited (if any) Hamiltonian control is 
available. 

In order to establish whether pd can be stabilized in the 
presence oi Hq, we can first check whether Hq\'^) = A|^), 
in which case the state is invariant. If not, the first step 
is to see whether invariance may be enforced by applying 
suitable QL dissipation. Consider a QL decomposition 
of Hq = J2k-^k' ^^^'^ decompose H^ in matrix blocks 



1. Determine whether pd would be DQLS in the ab- 
sence of Hq , by applying Theorem [1] 

2. If Pd is DQLS for Hq = 0, a straightforward mod- 
ification of Theorem [2] above (see also Corollary 
[5] in Appendix IVIip proves that a generic choice 
of stabilizing noise operators, {Dk}, will make pd 
GAS under the combined evolution generated by 
C{Hl„{D',}[J{Dk}). 

3. If Pd is not DQLS for Hq = and some QL Hamil- 
tonian control is available, we still invoke Theorem 
[5] to determine whether a stabilizing QL Hamilto- 
nian He (and possibly stabilizing noise operators 
{Dk}) can be found by randomizing the free pa- 
rameters in the controlled generator. If pd is QLS 
and solutions are attainable with the available con- 
trols, with probability one a choice will be found, 
such that Pd is GAS under the combined generator 
CiHl,+He,{D',}lJ{Dk}). 

4. If Pd is not DQLS and no Hamiltonian control is 
available, we need to check whether Hq satisfies 
the necessary conditions of Proposition[3] (or Corol- 
lary [Sj. If so, we invoke Theorem [2] to determine 
whether stabilizing noise operators {Dk} can be 
found by randomizing the free parameters in the 
controlled generator. If the drift Hamiltonian Hq 
does not not obey Proposition [3] and He = 0, pd 
cannot be made GAS. 



according to H = Tid < 



> "Hj-, that is: 



5. li Pd cannot be made GAS, conditional stabiliza- 
tion may be attempted, by invoking TheoremUand 
a suitable modification of the randomized stabiliza- 
tion algorithm. Specifically, in Step it is necessary 
to additionally ensure that the candidate subspace 
Ti' and its complement fi'-^ are invariant under the 
drift dynamics generated by Hq. By denoting with 
H', as before, the orthogonal projector onto H', this 
translates into requiring the additional compatibil- 
ity requirement [i7Q,H'] = 0. If pd can be made 
H'-DQLS with H' and "H'^ invariant, the algorithm 
succeeds with probability one. 



H" 



S,k 



P,k 



HpM 
HR,k 



For each fc, define: 



Dk 



1 Dp^k 




D 



P.k 



2iH 



p.k- 



(23) 



By Corollary 1 in [T2|, it follows that applying Dk as a 
noise operator for each neighborhood makes pd invariant 
for the global dynamics. Therefore, by recalling Eqs. ([5])- 
(|6]) and Corollary[l] we can find an equivalent representa- 
tion of the generator in standard form, say C{Hq, {D'f.}), 
with -ffol^) = and pd invariant under C. We may then 
proceed as follows: 



B. Drift Hamiltonian and dissipation 

In the most general situation, the system may be 
driven by a QL Markovian drift dynamics, specified by 
a generator £q{Hq,{L1}) with both Hq and some noise 
generators being non-vanishing. In this case, the first 
key property to verify is whether the necessary condi- 
tion for invariance of pd is obeyed by the dissipative 
drift, that is, whether for each neighborhood, we have 
Mk, Lk\^) = Afel*). If Lk\^) ^ Afel*) for some k, then 
Pd cannot be made GAS by the methods we described. 
While in practice one may expect that approximate sta- 
bilization be still meaningful and viable if the natural 
dissipation is sufficiently weak with respect to the avail- 
able controlled dissipation, establishing rigorous results 
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in this sense requires a separate analysis, which will be 
addressed elsewhere (see also [i^ for relevant distance 
bounds). 

Let us thus assume that Lk\'^) = Afe|^) for all k. Then 
by using Corollary [U we can write an equivalent genera- 
tor in standard form, with operators {Hq, L'f.} such that 
L'j,|^) = 0. By the invariancc requirement, these noise 
operators must satisfy Lemma [T] as well. In this way, we 
have effectively mapped the problem back the one just 
considered in Sec. IV A[ with only Hamiltonian drift. In 
case sufficient QL Hamiltonian control is available, we 
can remove the effect of Hq as above, and the problem is 
reduced to determine whether pd is DQLS, QLS, or H'- 
DQLS, as already discussed. Otherwise we may proceed 
as follows: 

1 . If the drift Hamiltonian Hq cannot be canceled by 
Hamiltonian control, we need to determine whether 
it destabilizes the desired state. If so, the action of 
Hq can be compensated in a way similar to the 
one given in Eq. ([251) . but, due to the non-zero 
dissipative drift {ijj}, we now need to choose 



D, 



1 Dp,k 




D 



P,k 



2iH 



P.k 



T'' T' 
^ S.k^P,k- 



The drift dynamics induced by TJg, {i'j,}, plus the 
QL operators {Dk}, can then be associated to a 
new generator £(iJ^, {L'i}) ^ C{H'^, {L'J U{^4) 
in standard form, with Hq\^) = and pd invariant 
under C. 

2. If Pd is DQLS for H'^ = 0, then by Corollary [i a 
generic choice of stabilizing operators, {Dk}-, will 
suffice to make pd GAS under the combined evolu- 
tion generated by /:(iJ^', {L'D Ul^fc})- 

3. If Pd is not DQLS for Hq = and some QL Hamil- 
tonian control is available, we can use Theorem [2] 
to determine whether a stabilizing QL Hamilto- 
nian He (and possibly stabilizing noise operators 
{Dk\) can be found by randomizing the free pa- 
rameters in the controlled generator. If pd is QLS 
and solutions are attainable with the available con- 
trols, with probability one a choice will be found, 
such that Pd is GAS under the combined generator 
C{H'^ + HeALl}[j{Dk}). 

4. If the state is not DQLS or Hq = and no Hamilto- 
nian control is available, we need to check whether 
Hq satisfies the necessary conditions of Proposition 
13] (or Corollary [3]) . If so, we are left with the prob- 
lem of finding stabilizing QL noise operators to be 
added to the internal ones, by randomizing the free 
parameters in the total controlled generator. If Hq 
has a destabilizing action on pd, the state is not 
QLS. 

5. If Pd is neither DQLS nor QLS, conditional stabi- 
lization may be tested, in a way similar to what de- 
scribed for Hamiltonian drift only. In Step of the 



stabilization algorithm, we now need to to ensure 
that both the candidate subspace H' and H'-^ are 
invariant under the full drift dynamics generated by 
Cq. If Pd can be made H'-DQLS with H' and "H'^ 
invariant, the algorithm succeeds with probability 



VI. ILLUSTRATIVE EXAMPLES 



Consider a quantum register consisting of n qubits, and 
let {(Ti \i = x,y, z} denote single-qubit Pauli matrices [l|, 
with the corresponding multi-qubit Pauli operators given 



by CT, 



(a) 



1 1 ® I, and at acting on the 



a-th qubit. Since the form of possible drift dynamics is 
highly system-dependent, we focus here on the drift-less 
setting, as relevant in particular to trapped-ion open- 
system quantum simulators [J| and cold atomic gases [1^ . 



A. DQLS states 

As mentioned earlier in the text, the DQLS states in- 
clude important classes of pure multipartite entangled 
states, most notably, all stabilizer and graph states, rel- 
ative to the natural choice of neighborhoods associated 
with connected nodes on the graph [Tlj, [3g ■ While char- 
acterizing the class of DQLS states for a pre- determined 
neighboorhood structure remains an interesting open 
problem, one can show that non-trivially entangled pure 
states which are provably not graph states are in DQLS. 
For instance, as already noted in [3^, the DQLS set for 
n ~ 4 qubits contains the symmetric Dicke state with 
two excitations. 



L>f ^ = (|0011) + |0101) + |1001) + |0110) + |1010) + |1100))/V6, 

which may be used as an optimal quantum resources in 
1 — > 3 telecloning protocols and has been experimen- 
tally characterized using polarization-entangled photons 
|39j. Imagine that the four qubits are arranged in a 

(2) 

line. The state D\ may be projected into both the 
three-qubit GHZ and the W class via single-qubit mea- 
surements. While the resulting GHZ3 and W3 states 
are graph states (hence DQLS) relative to the natu- 
ral ('star'-connected) neighboorhood choice A/i = {1,2}, 
JV2 = {2, 1, 3}, A/a = {3, 2}, this is no longer the case if 
a nearest-neighboorhood coupling constraint is imposed, 
in which case A^i = {1, 2}, 7V2 = {2, 3}. 

In what follows, we explicitly address dissipative sta- 
bilization of n-qubit GHZ and W states, which play a 
paradigmatic role in QIP [S^ ■ 
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B. Greenberger-Horne-Zeilinger cat states 

1. GHZ states are not DQLS 

A representative of the GHZ class is the state pghz 
|*)(vl/|, where 



I*) = |*GHZ> 

= (|0)®" + |1)®")/V2 EE (|000 . . . 0) + |111 . . . 1 



As shown in [36|, GHZ„ states are generally not DQLS, 
except in cases where the QL constraint becomes effec- 
tively trivial, as illustrated in the above three-qubit ex- 
ample. In fact, any reduced state on any (nontrivial) 
neighborhood is an equiprobable mixture of |000...0) 
and |111 . . . 1). It is then immediate to see that 

-Ho = span{|000 ... 0), |111 . . . 1)} = f]suMPM^ <^ l^J, 



with dim('Ho) = do = 2. From the point of view of MPS 
theory, the GHZ state is known to be non-injective, and 
corresponding to the two-fold degenerate ground state of 
a frustration- free (gapped) Ising parent Hamiltonian |34| . 



2, 



If we allow one neighborhood to include n^ = n/2 
qubits, we can always construct a QL Hamiltonian such 
that 

iJc(|000...0) = (|1...10...0) - I0...01...1 

iJc(|lll . . . 1) = (-|1 . . . 10 . . . 0) + |0 . . . 01 . . . 1 

with the vectors in the r.h.s. containing precisely n/2 
zeroes and n/2 ones, which clearly satisfies the require- 
ment. This may be achieved by considering a neighbor- 
hood structure consisting of two disjoint sets S£^i^2, each 
including half of the qubits if n is even [or (n+l)/2, (n — 
l)/2 in the odd-n case], and QL Hamiltonian components 
given by 



Hc = {-iY-'{l[-. 



(a) 



2. GHZ states can he QLS: conditions on the neighborhood 

size 

Let us explore under what conditions adding Hamilto- 
nian control can render pcHZ GAS. Given Corollary |3l 
we seek a QL Hamiltonian He = X^fc ^k such that: 

i/c|*GZH> =0, 

i/c|$i) i no, |$i) = (|ooo ... 0) - |iii . . . i))/V2. 

In order for this to happen, He must equivalently obey 
the condition 

iJc|000...0) = -Hc\lll...l). 

Since each component Hk is QL, it acts non-trivially on 
at most a number nk of "symbols" (that is, or 1) in the 
factorized states |0)®", |1)®", respectively. Let n be even 
and assume that nk < n/2. Then it follows that 

i7c|000...0) e spanUxi,...,Xn),Xj e {0,l},'^Xj < n/2\, 

3 

i7c|lll---l> e span||a;i,...,a;„),Xj G {0, l},^2:j > n/2\. 

3 

Thus, the two vectors must be orthogonal, since they 
belong to subspaces spanned by two orthogonal sets of 
vectors. This means that a Hamiltonian satisfying the 
requirements of Corollary |3] does not exists. In other 
words, we need He to be able to flip at least n/2 qubits 
in the product basis states |0)®", |1)®". 



ae5f 

If there exists two neighborhoods Nki , compatible with 
the QL constraints, such that Si C Nki, then, with a 
proper choice of the noise operators, a GHZ state may 
be rendered GAS. Otherwise, the above argument proves 
the following: 

Proposition 5 Assume that no neighborhood exists, that 
contains n/2 qubits if n is even, or (n + l)/2 if n is odd. 
Then pcHZ is not QLS. 

While the randomized QLS algorithm may be applied 
if a large neighborhood exists, stabilizing generators may 
be constructed by inspection in simple low-dimensional 
situations. Let us reconsider, for example, n — 3 qubits 
on a line, under the QL constraint of two-body nearest- 
neighbor couplings. We may then choose the two disjoint 
sets Si ~ {1},S2 ~ {2,3}, and implement 



Hr.=<J 



-Ji) 



.(2) 



r(3) 



i?l=/®(|00)(01| + |ll)(10|), 

D2 = /®(|00)(01|+i|ll)(10|). 

We stress that the phase factor appearing in the Lindblad 
operator D2 is not coincidental: in fact, the more sym- 
metric choice D^ = I® (|00)(01| + |11)(10|) would leave 
the — 1-eigenspace of the operator {^^^i 2 3 "^^^ invariant 
for the Liouvillian generator associated to Hc,Di,D'2. In 
constrast, pghz is the unique invariant state for the gen- 
erator associated to Hc,Di,D2, as required for GAS. 

This simple example is sufficient to establish that 
DQLS C QLS, namely that there exist pure entangled 
states that are not stabilizablc by dissipation alone but 
can indeed be made GAS by the addition of a suitable 
QL Hamiltonian. Furthermore, our argument for render- 
ing GHZ states GAS with QL resources is general, and, 
in the light of Proposition [51 it requires neighborhoods 
of the minimum possible size. Since the latter scales lin- 
early with the number of qubits n, the order of interac- 
tion in the required control Hamiltonians likewise grows 
extensively with n, preventing QLS to be achievable in a 
scalable fashion. 
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3. GHZ states are conditionally DQLS with a scalable 
strategy 

Consider again a three-qubit register, with locality no- 
tion associated to A/i = {1, 2}, and N2 = {2, 3}. In order 
to show that GHZ states can be H'-DQES for a properly 
chosen "H', note that |^)ghz is a -|-l-eigenvector of af^ , 
whilst |$i) = (|000) - |lll))/v^ is a -1-eigenvector of 
the same observable. In the notation of Section llVi we 
have T-Lyj = span{|<I>i)}. We can thus try to choose T-L' 
as the -l-1-eigenspace of crf^, denoted by 'Htxx- ^^ noise 
operators, pick D12 ~ D (E) I, and D23 ~ I ® D, with 



D 



10 





10 



<y. + I 



-, (24) 



in the standard Pauli formalism. It is easy to show 
that they both commute with cr®^ , and leave pcHZ in- 
variant. Finally, by running the DID algorithm for 
GAS verification described in Appendix lVIII starting from 
"Ho = span{|000), |111)} as the target subspace, one finds 
that it runs to completion in 4 steps and hence Ho is 
GAS. Hence, by Theorem [3l we conclude that pghz is 

^xa;x-DQLS. 

Does the same strategy still work for a generic number 
of qubits? More precisely, consider n qubits arranged on 
a linear graph (equivalently, an open spin chain), with 
neighborhoods associated to all nearest-neighbor pairs, 
namely A/fc = {fc, fc-|-l}fc=i^..._„_i. Implement a dissipator 
for each neighborhood of the form Dk ~ D ® Ij^ , where 
D has the structure given in (P^ . Consider, as in the 
above three-qubit case, any initial state with support on 
H', the -|-l-cigenspace of the operator af". Explicitly 
checking GAS of the n-qubit subspace Hq by resorting 
to the DID algorithm becomes impractical as n grows. 
Luckily, Theorem [3] makes this step unnecessary. In fact, 
note that: 

1. This choice of Dk makes 2)(H^ ) GAS on any 
neighborhood, so that S)('Ho) is GAS; 

2. The Dk commute with cr®" and leave %' invariant; 

3. The other state in "Ho, |$i) does not belong to T-L' . 

We can thus directly conclude that the GHZ„ state is 
"H'-DQLS, with the same choice of two-body dissipators 
for any n and the dimension of the conditional prepara- 
tion subspace being equal to d/2 = 2"~^. In practice, 
initialization to any state in T-L' by either coherent or 
incoherent means will suffice. For instance, two possi- 
ble scalable strategies are initialization into the product 
state I -I-)®" by application of a collective Hadamard gate 
7f®" to |0)®" (in case the latter is a natural starting 
point) or by a projective measurement of the collective 

(one-body) spin observable Sx = "^a"^^ post-selected 
on the outcome corresponding to the maximal eigenvalue 
(highest-weight) state. 



W states 



1. W states are not DQLS 



A representative of the W class is the state pvv = 
I *)(*!, with 

I*) = |*w) = (|100...0) + |010...0)+...+|000...1))/V^. 

This state is a permutation-invariant superposition of 
all computational basis states with a single 1, thus also 
a symmetric Dicke state one-excitation. The reduced 
states on any non-trivial neighborhood are statistical 
mixtures of |000 ... 0) and a smaller W state, say I^vvOj 
whose dimension is determined by the neighborhood. 
Accordingly, as established in [3g|, pw is generally not 
DQLS; since 

•Ho = span{|000. . .0), |*w)} = {~\^Vi]i]i{pM^ ® IatJ, 



with dim('Ho) — d^ — 2. Like the GHZ state, the W state 
is also known to be non-injective, however a frustration- 
free gap ped parent Hamiltonian does not exist in this 
case [3411. 



2. W states can be QLS: two-body interactions 

We here show that W states can be made QLS in 
the presence of QL constraints that prevent DQLS to 
be achievable. Consider, in particular, a n-qubit regis- 
ter with neighborhoods associated to arbitrary subsystem 
pairs, that is, Afjk = {j, fc}, j,k = 1, . . . ,n^ j ^ k. Since 
Ho — span{|\I'w), 10*^")}, we are in a situation where 
Corollary [3] applies. 

In order to construct a control Hamiltonian that sat- 
isfies the requirements of Corollary [31 it is useful to note 
that the target state may be rewritten as 



*w) = 



n-2 



|0)|*w„_.)|0) 



:(|1)|0)' 



K"-2) 



-|0)|0)®("-2)|1)). 



|0) 



It is then easy to show that, for example, the following 
two-body Hamiltonian satisfies the desired conditions: 



n-l 



H, - ^DPo - Po^i") , ^0 = ^ ^i"^ -in- 4)/('^) . 

a— 2 

This follows from the fact that Pq is an operator on 
0a=2 ^a which obeys 

Pol^w^^^) = 0, Po|0)«^("-') = |o)®("-2). 
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While we derived this particular Hamiltonian guided by 
simple symmetry considerations, to our scope it suffices 
to verify by direct computation that 



H,\^w) = 0, H,\0)^ 



I10...0)- I0...01), 



as needed for Corollary [31 In order to conclude that 
the above Hamiltonian indeed makes JV&vk) QLS, we also 
need to exhibit an explicit choice of dissipators Dk- To 
this aim, consider in each two-qubit neighborhood a lad- 
der dissipator of the form 



D 



110 

1 

1 





I(E)(J++(J+<E)I. 



This choice makes the set of states with support on H'^ 
GAS on each neighborhood, and hence Hq is stabilized. 
Explicit numerical calculation demonstrates that the 
above combination of Hamiltonian and dissipative con- 
trol makes \'^w) GAS, and hence QLS, at least for 
n = 3,4. While we have no general formal proof, we 
expect that the same protocol will work for arbitrary n. 
This shows how W states can be made QLS by allowing 
for arbitrary two-body interactions. Whether scalable 
stabilization protocols may be constructed in other rele- 
vant QL scenarios, for instance solely involving nearest- 
neighbors two-body interactions, remains open to further 
investigation. 



3. W states are conditionally DQLS with a scalable strategy 

Interestingly, W states can be conditionally stabilized 
by employing the systematic approach of Section IIVBI 

The key step is to show that the strict inclusion condi- 
tion given in Eq. (|14p holds. As we already noticed, 
with respect to any neighboorhood topology we have 
Ho = spanjl^vy), |0)®".} Let dk denote the number of 
subsystems in the neighboorhood Mk and, as before, let 
I'^Wa ) be a W state on dk qubits. We thus have to check, 
if for every k: 

supp(Tr^J|0)(0|®")) = span{|0)®'^n 






which is clearly true. Consider then 



■ATfc 



n'=nQf]nk = ■Hep|span{|0)®'^'=}«)-H 

k k 

= Hespan{|0)®"}. 



Physically, we may think of H' as the subspace of states 
orthogonal to the vacuum, thus initialization in "H' may 
be achieved in principle by any (coherent or incoherent) 
means that creates at least one "excitation" . 



"0 





... 








1 ••• 







1 ■•. 



In each neighborhood, consider now |0)®'''', |^^^) in 
this order and complete it to an orthonormal basis for 
T~l-j\fi^. With respect to this basis, define the following 
ladder-type Lindblad operator: 



Dm.^ 



The corresponding dissipative process may be thought as 
coohng the system to the reduced W state IvP^), while 
leaving the (QL) ground state |0)'^'^'= invariant. Then 
construct the overall dissipators as Dk = Dj^^ ® Ij^ . 
This choice ensures that: 

1. %' is invariant for each Dk\ 

2. SCHo) is GAS. 

Corollary [6] may then be invoked to establish that \^w) 
is "H'-DQLS for the choice of T-L' and {Dk} we made. 
Remarkably, the proposed protocol is both scalable and 
portable, in the sense it works for an arbitrary number of 
qubits and arbitrary QL notions. 



VII. CONCLUSION AND OUTLOOK 

We have provided a system-theoretic analysis of dif- 
ferent scenarios and strategies for designing Markovian 
evolutions that have a desired pure entangled state as 
their unique stable steady state, using time-independent 
control parameters and subject to realistic locality con- 
straints. In particular, we have shown how target states 
that are not stabilizable under purely dissipative QL con- 
trol, as previously considered in [36[, may be dissipa- 
tivcly prepared upon restricting the allowed set of ini- 
tializations or by allowing a combination of Hamiltonian 
and dissipative control. We have further addressed the 
role of Hamiltonian and/or dissipative drift dynamics, 
as arising from possible always-on coherent interactions 
and/or couplings of the target system to an uncontrol- 
lable Markovian environment. Constructive algorithms 
for synthesizing effective choices of stabilizing Markovian 
generators have been presented, suitable in principle for 
open-loop control implementations based on a switch- 
ing output-feedback law along the lines described in [3a] . 
In particular, if the Markovian semigroup is obtained as 
an average over the trajectories of a stochastic master 
equation, we recall that convergence of the semigroup en- 
tails convergence of the underlying stochastic dynamics 
in probability [44 1. 

While our results substantially expand the theoreti- 
cal framework and toolbox for QL dissipative entangle- 
ment engineering in QIP, a number of open questions 
and further directions for exploration exist. For time- 
independent Markovian dynamics as considered thus far. 
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issues of efficiency and robustness are especially impor- 
tant from a practical standpoint: (i) On the one hand, 
once a target state of interest is found to be stabilizable, 
it is desirable to characterize the speed with which asymp- 
totic convergence is attained depending on the system 
size and the elapsed (in reality always finite) stabiliza- 
tion time. While the analysis carried out in [ij] indicates 
that the relevant Liouvillian spectral gap scales favorably 
with n at least for injective MPSs, additional work is 
needed for more general classes of states as well as for 
optimally "tuning" the convergence speed as a function 
or the available control parameters, in the spirit of |35l |. 
(ii) On the other hand, assessing how sensitive steady 
states are with respect to deviations of the actual control 
parameters with respect to the intended ones, and/or to 
additional (static or time-dependent, possibly quantita- 
tively unspecified) perturbations is worth being pursued 
both in terms of analytical bounds and numerical explo- 
ration in specific QIP settings. General results recently 
established in [43, |4a] may prove useful in that respect. 
Partly related to the above, given that finite evolution 
time and errors will typically prevent a pure steady state 
to be exactly achieved, addressing the general problem of 
dissipative QL mixed-state stabilization is also a natural 
important next step. Additional physical motivation to 
consider target mixed entangled states is provided by the 
possibility to characterize and engineer non-equilibrium 
critical behavior in noise-driven many-body systems, so- 



called dissipative phase transitions, as recently investi- 
gated e.g. in [40|-|43J and already experimentally explored 
in [5| in small dimension. 

Finally, a yet different direction is provided by the in- 
vestigation of dissipative entanglement engineering in dif- 
ferent classes of open quantum dynamical models: from 
this point of view, both switched QL Lindblad dynam- 
ics generalizing on the work of [6[, and discrete-time 
Kraus-map engineering [43| are especially promising in 
the near term. Developing a stabilization framework for 
continuous-time non-Markovian dynamics is also an im- 
portant challenge down the line, as relevant to both un- 
derstanding the role of non-Markovianity as a resource 
[48J and to open-system quantum simulators |^, i4^ ,5^ . 
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Lemma 4 Let f{x) : 



pM 



function in x ^ R and let S = {x E M | f{x) 
Then /i(5) = 0, where fi is the Lehesgue measure in 



he a (non-zero) analytic 
0}. 

B-R" 



Define an m x n matrix X = [fjk{x)\, with fjk : "^^ — > 
C, such that the real and imaginary parts ?ii{fjk), '^(fjk) 
are (real)-analytic, and let r^ = maxj-gcx rank(X). No- 
tice that rank(Ar) S {0, . . . ,min{n,?7i}} for all x € C^ , 
and the maximum is attained in M.^ . We thus have the 
following: 



Lemma 5 The set X 

such that fi{X) = 0. 



{x e R''^ I rank(X) < r^} 



IS 



Proof: Let x £ R-'^ be such that X = [/jfe(a;)] has 
rank(Ar) = r,„. Then there exists a square r^ x r^ sub- 
matrix Xm of X that has rank r^ [33 ■ Notice that 
det{Xm) is a (complex) polynomial function of the el- 
ements of X, fjk, and hence of ^{fjk),'^{fjk)- Thus, 
5ft(det(X„i)) and 3(det(Ar„i)) are (real)-analytic func- 
tions of a: G R", and by virtue of the previous Lemma 
the set: 



X = 



{x e R^ I det(X™) = 0} = {xe 

= 3(det(X„))} 

{x e R^ I rank(X) < r„} 



oK 



5R(det(X„0) 



is the intersection of two zero-measure set. Thus, it itself 
satisfies fJ.{X) = 0. D 

: Dissipation-induced decomposition and randomized 
stabilization 

A general approach to decide the stability of a sub- 
space, or more precisely, of the set of states with support 
on a subspace, has been developed in [33 ■ We recall 
here some basic ideas and results, expanding the (par- 
tially randomized) strategy for Hamiltonian control de- 
sign proposed there to the synthesis of both coherent and 
dissipative controls. 

Let Hs be a proper subspace oiH and 'D{Hs) C J)('H) 
the set of states with support on Hs- Then 'D{Hs) is 
GAS for the QDS dynamics if and only if a Hilbert space 
decomposition in orthogonal subspaces, of the form 



n^Hs^n 



(1) 



i-H, 



(2) 



i-H, 



(?) 



(25) 
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can be obtained as the output of a constructive algorithm 
for GAS verification [33] (see also next section). Such 
decomposition is called the Dissipation-Induced Decom- 
position (DID). Each of the subspaces H^ in the direct 
sum is referred to as a basin. 

Partitioning each matrix associated to the noise op- 
erators Dk in blocks according to the DID results in 
the following standard structure, where the upper block- 
diagonal blocks establish the dissipation-induced, cas- 



cade connections between the different basins H 



Dk^ 



Similarly, for the control Hamiltonian we get: 



(') . 



' Ds 


b^^^ ••• 1 





4^^ 4^ ••• 




i?w < bf ■■. 







(26) 



Hr 



Hs 



H 



WF 



H 



(0) 







H. 



W 



(27) 



By construction, the bp blocks are either zero or full 
rank. The fact that the first column of blocks has only 
Ds 7^ is a necessary condition for the invariancc of 
Ti^Hs)- It follows that bf) ^ 0, otherwise ©(Hs) cannot 
be GAS. 

Consider now the target pure state pd = |^)('I'|, its 
support T-Ld = Tis, H^ — ^^' ^^"^ ^ = ^s ® ^R- ^^~ 
sumc that we are free to design H, {Dk} in parametric 
form, as given in Eq. ([TT]) . The DID provides a key tool 
for proving Theorem [2l which we reproduce below: 

Theorem 2 // there exists a choice of ajk, l3jk G I = 
[—7,7] that makes pd QLS, then almost any choice of 
ajk,Pjk G T that makes pd invariant, makes it QLS as 
well. 

Proof: We begin by observing that if a choice of con- 
trols Hc,{Dk] make D('Hs) GAS, then for any A € M 
the rescaled controls AiJc, {'v/ADfe} also make it GAS, 
since the total Liouvillian generator is linear in He and 
quadratic in Dk- Assume then that a certain choice of 
ctjk,Pjk makes 'S{'Hs) invariant: one can always rescale 
all ajk,Pjk by A small enough so that they are all in 
[—7, 7] and still obtain a viable solution. Hence, it is not 
restrictive to look for stabilizing parameters in a bounded 
interval. 

Next, assume that ajk,Pjk rnake 'D{'Hs) invariant. 
Since, as shown in Section IIIBl it is not restrictive to 
assume that He, {Dk} are in standard form, we may take 
Hs C kcT{Dk) for all k. The invariancc conditions then 
translate in linear constraints on ajk,Pjk, namely: 



^ajfcCTjfcl*) 



0, 



'0jkajk\^) ^ 0, (28) 



Since by hypothesis pd is GAS, there must exist at least 
a solution. If the solution is unique, the statement is 
trivial. If the solution is not unique, by linearity there is 
a subspace of solutions. If so, we may re-parametrize the 
admissible solutions in the free-parameters Ajk,(3jk € K- 
We now focus on a DID that makes pd GAS, and that 
will be proved to be generic. By construction, the ma- 
trix block decomposition of the matrix representation of 
Hc,{Dk} must be such that for each iteration, indexed 
by j, cither 



f)U} 



D 



(3) 
PA 



D 



U) 

PM 



(29) 



has maximum rank r^^' — max{dim('Hg? ),dim('Hj? )} 
(see Step 2, 3. a and 3.b), or 



^P^=^Hf^\T.D'SlD^i 



(30) 



has full rank r^-'' (defined as above, see step 3.c of the 
algorithm). At the j'-th iteration, dim(Hj? ) is fixed, 
but we can choose a set of parameters Ujk , 13 jk that max- 
imizes dim(?^y^). Given Eq. (I29l)-(|30l), this is equivalent 
to maximize the rank of either V'^^^ or Cp at each it- 
eration. The elements of 2?(j) and Cp are (complex) 
polynomial functions of the real parameters ajk,Pjki and 
hence Lemma [5] ensures that the parameter set corre- 
sponding to maximal rank of the corresponding matrices 
has measure 1. This implies, in turn, that the same DID 
is constructed with probability one when the parameters 
are randomly chosen as above, establishing that the tar- 
get can be made GAS with probability one, as claimed. 
D 

Along similar lines, we may show that the set of pa- 
rameters that makes pd invariant has measure zero. In 
fact, such parameter choices correspond to the proper lin- 
ear hyperplane determined by Eq. (|28p . The same holds 
for the DQLS property: the proof is identical to the one 
above, just consider the ajk fixed. Thus, if there exists 
a choice of I3jk that makes pd DQLS, then almost any 
choice of (ijk, such that pd is invariant, makes it DQLS 
as well. 

Corollary 8 Assume He = X]fc ^k to be fixed, that is, 
ajk are given for all j, k. If there exists a choice of Pjk G I 
that makes pd DQLS, then almost any choice of jijk G I, 
such that Pd is invariant, makes it DQLS as well. 

As discussed in the main text, the randomized ap- 
proach carries over to conditional stabilization also, pro- 
vided that care is taken in ensuring that no quadratic 
constraints appear. While imposing invariancc of both 
Ti' and its complement H'^ is a mathematically natural 
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restriction to circumvent this issue, it need not always 
be easy to ensure in physical systems, especially in the 
presence of drift dynamics. Thus, an interesting open 
question is to determine whether a randomized approach 
may be devised for conditional DQLS in full generality. 

: DID algorithm 

In order to make our presentation self-contained, we 
reproduce here the algorithm for the construction of the 
DID [35|. The initial (target) subspace Hd ^ "Hs (or, 
better, ^CHs)), is GAS if and only if the algorithm runs 
to completion. 

The inputs are a QDS generator C specified by the 
operators H, {Lk} and an invariant subspace Us or the 
evolution, that is, such that e^*S)(Hs) C 1){'Hs)- Check- 
ing whether the above property holds can be easily done 
using Proposition [TJ 

Algorithm for GAS verification and DID construction 



Let Hs be invariant. Call H 



(0) 



Hr, H 



(0) 



ns, 



choose an orthonormal basis for the subspaces and write 
the matrices with respect to that basis. Rename the ma- 



r(0) 



trix blocks as follows: Hg'' := Hs, Hp' := Hp, H^ 



r(o) ., 



H 



= Hp, 



(0) ._ 



fl' ^SM ■— ^s,k, i^p^k ■— ^P,k, and L^^. 

For j > 0, consider the following iterative procedure: 

(1) 
1. Compute the matrix blocks Lp)^ according to the 



decomposition "H*^^ = WJ' S) H 






b. If n 



(i+i) 



R 



,0+1) 



^ {0}, but H 



a-Hi) 

R 



£ n'l"^ 



de- 



fine Tiqi as the orthogonal complement of 



H^i'+^^ in H^i\ that is, n^i+^^ 



n 



0) 



n 



0+1) 



c. un 



0+1) 



n^£ (that is, L^^l 



OVfc), define 



3^ := -tH'd^ 



~<j) 



-Ij:^ 



0)trO) 
k^R.k- 



If 0^> ^ 0, re-define U^^ 



0+1) 



If -H 



0+1) 



{0}, define U. 



0+1) 



ker(£ 



(j)^ 
P )■ 



n 



0) 



and the iterative procedure is successfully 
completed. Otherwise define H 

n^i^ e n 



0+1) 



■-fl 



0+1) 

R 



-lie 



0) 



,0) 



0, then T-Lp is invariant and T-Ls 



cannot be GAS. Exit the algorithm. 



4. Define -^1,^'+^^ := uf (SU^r^^ 



/O- 



-1) 



To construct a basis 



for T-Lg , append to the already defined basis for 

Hg an orthonormal basis for H^ ' . 
5. Increment the counter j and go back to step 1. 



2. Define V. 



0+1) 

R 



HfekcrLl^^,. 



3. Consider the following three sub-cases: 

a. If ■H^'+^^ = {0}, define 7^^+^^ := H^^ The 
iterative procedure is successfully completed. 



The algorithm ends in a finite number of steps, since at 
every iteration it either stops or the dimension oiUp is 
reduced by at least one. 



